Causal signal transmission by quantum fields. 
VII: Electromagnetic interaction of distinguishable devices and the optical equivalence 

theorem. 



O 

(N 

Oh' 
ON 



9 L,: 

c ■ 

=5 

cr 



> 
oo 
On 
O 

ON 

o 



X 



L.I. Plimak 1 and S.T. Stenholm 1 ' 2 > 3 

1 Institut fur Quantenphysik, Universitat Vim, 89069 Ulm, Germany. 
2 Physics Department, Royal Institute of Technology, KTH, Stockholm, Sweden. 
3 Laboratory of Computational Engineering, HUT, Espoo, Finland. 
(Dated: September 20, 2011) 

Generalised phase-space techniques for electromagnetic interactions beyond the rotating wave 
approximation [L.P. and S.S., arXive:1104.3825 (2011)] is applied to interactions of distinguishable 
devices. The paper is built around the concept of "doing quantum electrodynamics while thinking 
classically," which is a generalisation of Sudarshan's renowned optical equivalence theorem [E.C.G. 
Sudarshan, Phys. Rev. Lett. 10, 277 (1963)] to the interacting electromagnetic field. This concept 
allows one to reduce inherently quantum problems to semiclasical considerations. 

PACS numbers: XXZ 



I. INTRODUCTION 

This article concludes the subseries of this series of 
papers, concerned with quantum electrodynamics under 
macroscopic (often termed mesoscopic, cf. endnote 21 
in paper [l|) conditions. In paper Q, we applied re- 
sponse transformations and the causal Wick theo- 
rem [l| to the standard perturbative approach of quan- 
tum field theory Q. An astonishing feature of formulae 
thus found is that they lack Planck's constant. Such for- 
mulae survive the classical limit h — > unchanged, and 
must therefore equally make sense in classical stochas- 
tic electrodynamics. This quantum-classical correspon- 
dence was the subject of paper It was demonstrated 
that formulae of Ref. Q are naturally written in phase- 
space terms. We introduced the concept of conditional 
functional quasi-probability distribution, or conditional 
P- functional, where "conditional" means dependent on 
external c-number sources. The conditional P-functional 
generalises the conventional P-function [11| in two ways: 
to general nonlinear non-Markovian quantum systems, 
and to response properties of quantum systems. It also 
establishes a natural relation to classical stochastic elec- 
trodynamics. If a P-functional is nonnegative, it may 
be interpreted as a functional probability distribution, 
thus mapping a quantum system exactly on a classi- 
cal stochastic system (recall that corresponding formulae 
lack Planck's constant). 

The subject of Refs. @, was a solitary quantum de- 
vice. Formally, the latter is just a placeholder for a quan- 
tum model of matter. There are no restrictions on this 
model, so that results of Refs. @, @] apply to any case of 
electromagnetic interactions, including relativistic quan- 
tum electrodynamics. The macroscopic (mesoscopic) ap- 
proximation is introduced in this paper, by regarding the 
matter as a collection of distinguishable devices. Dis- 
tinguishability is an approximation: after all, electrons 
in all devices are identical. Under macroscopic condi- 
tions, such effects are negligible for all practical purposes. 



The approximation of there being distinguishable devices 
may be regarded a definition of macroscopic conditions 
in quantum mechanics. 

The notion this paper is built around is "doing quan- 
tum electrodynamics while thinking classically." Any re- 
lation for P-functionals is identical with some relation for 
probability distributions, and may be obtained as formal 
generalisation of the latter. One may therefore derive 
quantum relations by "demoting" problems to classical 
mechanics, doing the classical theory, then "promoting" 
the results to quantum mechanics. This recipe must 
be applied with care, because there exist classical mod- 
els incompatible with quantum mechanics, such as, for 
instance, noiseless photodetector and noiseless coherent 
amplifier. The necessary reservations constitute an im- 
portant part of our analyses. 

As in our previous papers [l|, 0, 0] we distinguish 
the narrow-band and broad-band cases, which differ in 
whether the rotating wave approximation (RWA) is or is 
not made in dynamics. The broad-band case is most 
general, but is also very much disconnected from the 
quantum-optical paradigm where RWA is quite common. 
For instance, in Glauber-Kelley-Kleiner's photodetection 
theory [§-floj , the optical field is commonly treated under 
the RWA. However, the latter cannot be extended to the 
photocurrent and photovoltage. These are broad-band 
processes which are not subject to any kind of resonance 
approximation. 

The paper is structured as follows. In Sec. HH we sum- 
marise the key results of papers @, 0|- I n Sec. IIII1 we 
demonstrate that the concept of "doing quantum elec- 
trodynamics while thinking classically" naturally extends 
to interactions of distinguishable devices. As an illustra- 
tion, in Sec. IIV1 we apply this concept to the so-called 
cascaded systems [13]. We formulate a photodetection 
theory without the rotating wave approximation (Sec. 
IIV C[) and a generalisation of Sudarshan's optical equiv- 
alence theorem to interacting electromagnetic field (Sec. 
IIV E[) . In Sec. IIV Fl we apply "doing quantum electrody- 
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namics while thinking classically" to photodetection of 
the electromagnetic field in a quantum state preceded by 
coherent quantum amplification. This example allows us 
to illustrate the reservations one has to make when ap- 
plying this concept to quantum systems. In appendices, 
we concern ourselves with formal particulars neglected in 
the main body of the paper. 



This definition is Kubo's formula for a linear response 
function [13| : for more details see Ref. Q. The com- 
mutator in ([6]) is a c-number so that quantum averag- 
ing present in Kubo's formula could be omitted (i.e., re- 
sponse of a linear system does not depend on its state). 
The explicit expression for Gr follows by combining (j6]) 
with ©, see Ref. 0. 



II. SOLITARY DEVICE REVISITED 



Condensed notation 



A. The model 

We start our analyses from the broad-band case [l|, 
HI 0] ■ We employ a somewhat simplified version of the 
structural model of electromagnetic interaction used in 

Refs. @, a. For convenience of the reader, we copy here 

the key definitions from Ref. [3] ■ We consider a quantum 
device interacting with a collection of oscillator modes, 
with the Hamiltonian in the interaction picture being, 



N 



(1) 



The oscillators, represented by the standard creation and 
annihilation operators, 

[o K) o£,] = 8 kK >, k, k' = 1, • • • ,N. (2) 
are organised in a quantised field, 



A(x,t) 



K.— l 



He. 



(3) 



where Uk(x) are complex mode functions. Variable x 
comprises all field arguments except time. The field inter- 
acts with the device according to the nonresonant (broad- 
band) Hamiltonian, 



Hj(t) 



dx[A(x,t) + A c (x,t)]j(x,t). 



(4) 



The Hamiltonian -ffdev(^) and the current operator J(x, t) 
describe the device. They commute with all o K ,a[ and 
otherwise remain arbitrary. The Heisenberg density ma- 
trix factorises into the vacuum state of all oscillators and 
an arbitrary state of the device, 

p= |0)(0| (gipdev (5) 

The c-number external source A e (x,t) is added for for- 
mal purposes. The aforementioned simplification is the 
absence of an external c-number current in the interac- 
tion which is of no use in this paper. 

B. Retarded Green function of the field 

The electromagnetic field enters the theory through its 
retarded Green function 

G n (x, x', t-t')= l -6{t - t') [A(x, t), A(x', t')) . (6) 



To keep the bulk of formulae under the lid and make 
their structure more transparent, we make extensive use 
of condensed notation, 



fg = J dxdtf{x, t)g(x, t), (7) 

fKg = J dxdx'dtdt'f(x,t) 

x K(x,x',t-t')g(x',t'), (8) 

(Kf)(x,t)= [ dx'dt'K(x,x',t-t')f(x',t r ), (9) 



(fK)(x,t) = J dx'dt'g(x',t')K(x',x,t' -t), (10) 

where f(x, t) and g(x, t) are c-number or q-number func- 
tions, and K{x, x' , t—t 1 ) is a c-number kernel. The "prod- 
ucts" fg and fKg denote scalars, while Kg and fK 
functions (fields). 



D. 



Conditional time-normal averages and 
conditional P-functionals 



A formal solution to the problem of Sec. Ill Al is con- 
structed applying the standard perturbative techniques 
of quantum field theory 0. Of interest is however not 
this solution as such, but the fact that there exists a 
"language" in which it looks essentially classical. Refer- 
ring the reader for details to Refs. 0, 01 , here we only 
reiterate the key points. 

As was demonstrated in Ref. 0] , full quantum descrip- 
tion of an electroma gne tic device reduces to time-normal 
averages [3, 0, fsl- llCf of the quantum fields and cur- 
rents, conditional on the sources. For the broad-band 
field and current, one should use the amended definition 
of the time- normal ordering introduced in Ref. 4] . The 
conventional as well as amended definitions of the time- 
normal ordering are reiterated in appendix IA3I For a 
general discussion of the time-normal operator ordering 
see Ref. Q, section II, and references therein. 

The said averages enter the theory through their gen- 
erating functional, 

$(r),C\A e )=(T:ex.p(ir]A + i{j):), (11) 

where r)(x, t) and £(x, t) are auxiliary c-number functions. 
Equation (|TT|) implies condensed notation ([7]). The oper- 
ators A(x, t), J(x, t) are A(x, t), J(x, t) in the Heisenberg 
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picture. These operators are by construction dependent 
(conditional) on the c-number source in the Hamiltonian. 
The averaging in (fTT|) is over the initial (Heisenberg) state 
of the system 

(■■.)=Trp(...), (12) 

where the ellipsis stands for an arbitrary operator. 

Formulae relating functional (fTTj) to conventional 
quantum averages of the Heisenberg operators (termed 
response transformations) are summarised in paper Q, 
section IIIB. In the terminology of Refs. @, 0, H, 0], that 
functional contains full information on the quan- 
tum device is a consistency condition. For a summary 
of consistency conditions see Ref. 0], section HID. Some 
examples may be found in Sec. Ill Fl below. 

E. Reduction to quantum current 

In classical stochastic electrodynamics, a device may 
be seen as a random current J(x, t). This current radiates 



the random field, 



A(x, t)= dx'dt'G R (x, x', t - t')J(x', t'), (13) 



where G R (x, x' , t — t') is the retarded Green function 
(also known as transfer function, or linear response func- 
tion) characteristic of the linear media, or vacuum, in 
which the device is submerged. Linear response func- 
tions of a classical field and of the corresponding quan- 
tised field coincide [!, [l4[ , so that Gr in (TT3")) is in fact 
given by Kubo's quantum formula (JBJ). 

If the radiating current is stochastic, joint statistical 
averages of the field and current are given by the formula, 



A\X\ , t\ ) • • • A(x m , tfn ) J (x m _|-i , t m +i )•••(/ (x m -|- n , ^m+n ) 

= J dx^dt^ ■ ■ ■ dx' m dt' m G R (x 1 ,x' 1 ,t 1 - t[) ■ ■ ■ G R (x m ,x' m , t m - t' m ) 

X J (x'i , t]J • • • J (x' m , t' m )J (x m +i , f m _|_ij ■ ■ ■ J (x m ^- n , t m + n ) , (14) 

As was shown in Ref. Q, the corresponding quantum formula for joint time- normal averages of the field and current 
operators generated by functional (jlip is found replacing statistical averages by time-normal averages, 

(^l~'A{x\ , t x ) • • • ^A(^C m , tjy^^J (x m _|_i , ) ' ' ' 3 (^m+n j ^m+n ) •} 

= J dx\dt' 1 ■ ■ ■ dx' m dt' m G R (xi,x\,ti - t[) ■ ■ ■ G R (x m ,x' m ,t m - t' m ) 

x {T'-j(x'i , t'i) ■ ■ ■ j(x' m ,t' m )j(x m +i, t m +l) ' ' ' <J{x m +m Wn):). (15) 

I 



One may say that, under the time-normal ordering, Eq. 
(fT3"| applies directly to Heisenberg operators. 



F. Response and phase-space characterisation of 
solitary devices 

1. "Dressed" device 

We remind that, by definition, J{x,t) is the 
interaction-picture operator ( "bare" current). Its Heisen- 
berg counterpart ("dressed" current) is denoted as 
J(x,t); it is by construction dependent (conditional) on 
the external field A e (x,t) present in 

Equation (|15l) has an obvious implication: it suffices to 
calculate time-normal averages of the quantum current. 



Those of the field are recovered applying the classical 
radiation law (|13l) . Using Eq. (|15l) . for the functional 
(flTt we obtain, 

$(r),<;\A e ) = <P(0,( + vRr\A) = <P dcv (( + v R R \A c ). 

(16) 

The "dressed" device is completely characterised by the 
time-normal averages of the Heisenberg current operator, 
generated by the functional, 

*dcv(C|^c) = (T:exp(iCJ):) 

= Il{/ dJ(x,t)\p(j\A e )exp 

(17) 
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We also took this opportunity to introduce the phase- 
space characterisation of the "dressed" current by 
the conditional P- functional , or conditional functional 
quasiprobability distribution, p(J\ A c ) 7]. 



2. Consistency condition 

An alternatively way of defining functional (| 1 T[) is 
applying response transformation [3m| to the closed- 
time- loop (Schwinger-Perel-Keldysh [15T - ll7l |) averages of 
S{x,t) defined without the source. Namely, 



dov(C|a ) = [(Tcexp(iC+J+-K-J-))\ Ae=0 ]\ c 



where c.v. refers to the response substitution, 
a c (x,t) 



C±(x,t) 



±C {T) (x,t), 



(18) 



(19) 



with W standing for separation of the frequency-positive 
and negative parts of a function. Definition of the Tq- 
ordering is reiterated in appendix IA H and that of the said 
separation — in appendix IA 21 

Equivalence of definitions (ITT1) and (fT5)l is, in termi- 
nology of Refs. 011,0], a consistency condition. More 
precisely, the latter is expressed by the relation H| , 

$dov(Ch + A c ) = (T c exp (i( + J+ -iC-J-))\ cv . 

(20) 

showing that the auxiliary variable a c (x,t) and the ex- 
ternal source A c (x, t) occur in the theory as a sum. With 
a e (x, t) = we recover Eq. (|17|) . while with A c (x, t) = 
- Eq. (|18[) . A summary of consistency conditions 
may be found in Ref. Q, section HID. For details see 
Refs. 0,110 

Equation (f2~4l gives a fair idea of how the language of 
conditional time-normal averages is related to the con- 
ventional closed-time- loop formalism [l5l - [T7l |. In partic- 
ular, it makes it evident that the time-normal averages 
(fTTf indeed provide complete quantum characterisation 
of the device. 



3. "Bare" device 

Similar to the "dressed" current, the "bare" current 
may be characterised in two equivalent ways. One is in 
terms of the time-normal averages of the Heisenberg cur- 
rent operator in the presence of a given c-number source 
with the quantised electromagnetic field "switched off." 
These are conveniently accessed through their generating 
functional, 



$ dcv (C| A c ) = TrpdevT : exp (i(j') 



m / dJ(x,t)\ P 1 (j\A c )exp(iCJ), 



where J'(x, t) is defined as a Heisenberg operator for the 
Hamitonian, 

H\t) = H dcv (t) - J dxJ(x, t)A c {x, t), (22) 

and ^(jjAe) is the corresponding "bare" P-functional. 
Alternatively, 

$dcv(C|ac) = (Tcexp (*C+J+-iC-J-)>| c . v .. (23) 

Equation (|23p makes it evident that $ dcv and p l dcv are 
determined solely by the free current operator J{x,t). 
Equations (f2"Tj) and (f2"5| are particular cases of yet an- 
other instance of consistency condition, 

^ dev (C\a e + A e )=(T c exp(iC + J' + -iC-J'-))\ , 

C ' V (24) 

In fact, Eq. (|20l) reduces to (|2"4"|) if formally regarding the 
interaction with the quantised field as part of H dov (t). 



G. How to do quantum electrodynamics while 
thinking classically 

The main advantage of the quasiprobability distribu- 
tions is that they may be manipulated to a large extent 
as if they were classical probability distributions. For ex- 
ample, one may introduce a joint quasiprobability distri- 
bution of the quantum field and current by the formula, 



T: exp 

= n|y dJ(x,t)dA(x,t)\p(A, j\A c )exp(inA + i(j). 

(25) 

In classical stochastic electrodynamics, we would apply 
Eq. (TT3")) , resulting in, 



p(A, j\A e ) =p(j\A e )l[6(A(x,t) - (G R J)(x,t)), 

X,t 

(26) 

where we use notation ^ . Substituting this relation into 
Eq. dUD we find, 



r : exp (in A + : 
= []{/ ^,0}p(^|^c)exp(<J + l ? 7 G , R t /). (27) 



We use here notation ©. In view of Eq. pT|) . we have 
recovered Eq. (TT6")) . proving that Eq. (|2"6"|) is in fact a 
genuine quantum formula. 

Furthermore, the solution to the self-action, or elec- 
tromagnetic dressing, problem in terms of the quasidis- 
tributions reads, 



(21) 



P(J\A C ) =p L (j\A c + G R J), 



(28) 
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Jb 



b) 
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FIG. 1: Schematics of electromagnetic interaction of distin- 
guishable devices, a) Two solitary devices in response repre- 
sentation, b) The same devices interacting. 



where we use notation ^ . The classical content of this 
relation is obvious. Bare devices are characterised by 
statistics of the current conditional on the local field, 

Ai(x,t) = A c {x,t) 

+ ( dx'dt'G K (x,x' ,t-t')J(x' \t'). (29) 



Equation (|28|) states that, 

p(j\A c ) = P 1 (J\A 1 ) 



(30) 



As was shown in [7j , Eq. (f2"5f is equivalent to the operator 
dressing formula, 



$dcv(C|a ) = exp 
found in Ref. [2j|. 



i-^-G K -^)$ dev {(\a e ), (31) 



III. ELECTROMAGNETIC INTERACTION OF 
A PAIR OF QUANTUM DEVICES 

A. Statement of the problem 

The crucial step to physics is from a solitary device 
to a pair of interacting distinguishable devices. Assume 
that the quantum device consists of two components, A 
and B. Of interest to us is the connection among three 
quantum problems: those of solitary devices A and B 
(fig- Ut) ! and that of the composite device (fig. Hb). All 
three are governed by the Hamiltonian (JT|). The problem 



of composite device emerges by postulating, 

H dev (t) = H dcv A(t) + H dcvB (t), (32) 
J(x,t) = J A (x,t) + J B (x,t), (33) 

Pdcv = PdcvA ® pdcvB- (34) 

This implies factorisation of the matter subspace of the 
Hilbert space, so that H dovA (t) 7 J A (x,t) and pdcvA com- 
mute with H devB (t), J B (x,t) and p devB - Furthermore, 
the problem of device A is found specifying, 

H dev (t) = H devA (t), (35) 
J(x,t) = J A (x,t), (36) 

Pdcv — PdcvA- (37) 

while the problem of device B implies that, 

Hdev(t) = HdcvB(t), (38) 

J(x,t) = J B (x,t), (39) 

Pdcv = PdcvB- (40) 

Quantum averages defined under conditions (l32j) -(f34" |) . 
(jnSJ)— dSZI) and (JHEJ)— (SOJ) will be denoted, respectively, as 
(■■■),{■■■) A and 

Definitions of Sec. Ill Fl are extended to devices A and B 
by assigning subscripts A and B to functionals, averages 
and quasidistributions. So, 

$ dcvA,s(C|^c) = Tr/5 d cvA,i3T:exp (i(j' A , B ): 

= I1 1 / dJ A,B{x,t) [p\ B (J A , B \A C ) exp (<Ja,b), 
x.t ^ J > 

(41) 

where J' A B {x, t) are the Heisenberg currents with respect 
to the Hamiltonians, 



H(t) = H devA , B {t) - J dxJ A>B {x,t)A e {x,t). (42) 

Similarly, for dressed components, 

$dcvA,s(C|^c) =(T:exp(iC7):) 

= IT] / dJ A,B(x,t) \p A , B {J A , B \A e ) exp (i(J A , B ). 

(43) 



(44) 



Finally, eq. (|28|) turns into a pair of relations, 

p A (j A \A c ) = p\(j A \A c + GrJa), 
Pb{Jb\A c ) =p B (j B \A c + G R J B ). 



Assigning indices to phase-space variables in (|41l) - (|4"4")) 
is a matter of physical clarity rather than mathematical 
necessity. 
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Phase-space approach to interacting quantum 
devices 



It is not difficult to derive the formula relating $dev 
to QdcvA.B directly, using Eqs. (|3ip and the obvious re- 
lation, 



We use notation ©. We then dress the components of 
the device by taking notice of Eqs. (14"4"|) and (|4"6"|) . 



$Lv(C|^l) = $devA(CN$LvB(C|A 



(45) 



Such derivation is outlined in Sec. IIII Dl below. Our goal 
is, however, not to construct a formal solution to the in- 
teraction problem — which at this level of abstraction 
is more or less trivial — but to show that this solu- 
tion is consistent with "doing quantum electrodynam- 
ics while thinking classically." We therefore start from 
constructing the said solution in terms of the quasiprob- 
ability distributions, manipulating them as if they were 
classical probability distributions. In Sees IIII Cl and llTTDl 
we demonstrate that the formula thus found is equivalent 
to the exact q-number formula. 

The total random current characterising the composite 
device is a sum of two components, 



J(x,t) = Ja(x, t) + J B (x,t). 



(46) 



Device A is formally described by random current J a, 
which is characterised by two conditional quasiprobabil- 
ity distributions, p^(Ja|^4i) and pa(Ja\A b ) . For device 
B, we have current J B , characterised by p l B (j B \A\) and 
Pb(Jb\A c ) . The quasiprobability distributions are pair- 
wise connected by Eqs. (|4"4"|) , 

Solution to the interaction problem in terms of the 
"bare" distributions is trivial. The currents J a and J B 
are independent if considered conditional on the local 
field. The corresponding joint quasiprobability distribu- 
tion factorises, 



p 1 ( J A , J B I A x ) = pa {Ja \ Ai)p B (J B \Ay 



(47) 



Integrating over the redundant information we find the 
quasiprobability distribution of the full current, 

I?(J\Ai) = Jjjy dJA{x,t)\p l A{J A \A x )p B (J B \Ax), 

(48) 

where here and hereafter till the end of the paragraph 
J B (x,t) stands for the current difference, 



J B (x,t) = J(x,t) - JA{x,t). 



(49) 



Based on this relation it is also easy to construct a solu- 
tion to the interaction problem in terms of the "dressed" 
distributions. Applying the dressing formula ([25)1 to (|4"5)) 
we obtain, 

p(j\A c ) = \[^J dJ A (x,t)^p l A (jA\A c + G K j) 

x Pb (J b \A c + G r J). (50) 



(51) 



P 1 A (J A \A C + G R J) 

= p l A (J A \A c + G R [J A + Jb]) = pa(Ja\A cA ), 
Pb (J b \A c + G r J) 

= p l B (J B \A C + G R [J a + Jb] ) = Pb(Jb\A cB ) , 

where 

A cA (x,t) = A e (x,t) + / dx'dt'G R (x,x',t- t')J B (x,t'), 



A cB {x 1 t) = A e (x,t) + J dx'dt'G R (x,x',t-t')JA(x,t'). 

(52) 

Finally, 

p(j|A c )=J|i / dJA{x,t)\p A (JA\A c A)pB(JB\A cB ), 
x,t ^ J ' 

(53) 

The physical content of Eq. (|53l) is illustrated in fig. 
[Tp. It describes a random current J(x, t) which is a sum 
of two components, Ja{x,£) and J B (x,t) . Statistics of 
the components depend (are conditional) on the external 
fields A e A(x,t) and A cB (x,t). With these fields given, 
components of the current are statistically independent. 
Because of the electromagnetic interaction, the external 
fields themselves become stochastic. Each is a sum of 
the external field affecting the composite device plus ra- 
diation of the other component. From the point of view 
of the device, radiation of the components is part of the 
local microscopic field. From the point of view of each 
of the components, this radiation is part of the external 
field. Electromagnetic self-action of the components is 
included inp^B and is therefore excluded from Eq. (|53l) . 
In quantum mechanics, these considerations apply with 
replacement of "statistical" by "quasistatistical." The 
last statement remains subject to independent verifica- 
tion of Eq. ((53J in Sees UTTC] and HID] below. 

C. From quasiprobability distributions to 
generating functionals 

To assure a natural connection with quantum electro- 
dynamics we reformulate the theory of Sec. IIII 51 in terms 
of the generating functionals (for all definittions see Sec. 
EH). Substituting Eq. (gSJ into Eq. flST) we find, 

$Lv(C|^l) =n{ / dJ A (x,t)dJ B (x,t) 

x,t J 

x ^ a {J a \A^^ b {J b \Ai) exp [i((j A + Jb)]- (54) 

The integral factorises; using Eqs. (|4"Tj) we arrive at Eq. 
(05). 
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Finding the Hilbert-space counterpart of Eq. (|53l) takes 
a bit of ingenuity. Substituting Eq. (j53")l into ([XT]) yields, 



$dcv(ci^c) -n{ / 



dJ A (x, t)dJ B {x, t) 



x p A (j A \A c + G r Jb)pb(Jb\A c + GrJa) 

x exp [<(Ja + Jb)] ■ (55) 



Unlike Eq. (|54j) , here the integral does not factorise. 
However, we can write, 



Pa(Ja\A c + G R J B ) 



exp ^j^-G r Jb^pa(Ja\A c ), (56) 



and 



exp (<J B ) exp (j^- G R J B^j 



s s 

e-M> ( - i J^~ GR J^ ) '^i' K'^'. 1 



(57) 



Equation (|56p is an application of the functional shift 
operator, while (|57l) is just obvious. Relations similar 
to (I56p . (|57|1 may be written for the remaining two fac- 
tors in the integrand in (|55[) . Pulling the exponentiated 
differential operator out of the integral we find, 

(5 5 S 5 \ 
~ l w c GR rc' l sA c GR 6c) 

x Ii\ / dJ A (x,t)dJ B (x,t)\p A (J A \A e )p B (j B \A' e ) 

x,t ^ J J 

x exp (i(J A + i('J B ) | c , =Cijl , =Ae . (58) 

Introducing pairs of variables ((x, t), C( x , t) and 
A c (x, t), A' c (x, t) allows all differentiations to hit the right 
targets. The integral in (1551) is already factorised. Re- 
calling (|4"3")l we arrive at the relation sought, 



dcvS 



(59) 



D. Direct derivation of Eq. (|59 

Substituting (|35|> into (|3Tj) we obtain, 



*dov(C|a c ) = exp 



i— G R - 

5a e 5£ 



*dcvA(C|ae)$dcvs(C| 



(60) 



We now apply the relation [T3 |. 

= *(jf + Jrh {f)n{f %--,- ,61) 

where J r i(-), ^(O) -^(O are arbitrary functionals and 
f(x,t),f'(x,t) are auxiliary functional variables. Equa- 
tion (|6ip is a compact way of formulating general rules 
of product differentiation. It may be verified expanding 
•7"i(')j -^G); ^"bCO m functional Taylor series. Using it 
we rewrite (l60l) as, 



$dcv(C| 



exp 



5a 



Sa'J R UC *C 



X$ devA (C|ae)$ de vfl(C>e)L, , , ■ (62) 



Expanding the bilinear form in the exponent we have 

=exp (" i ^ GR l) exp (~^ GR l 7 
x exp ( - *£ Gi 4) cxp ( _ ^ Gr ^) 

The last two factors here "dress" the devices, 



(63) 



exp 

exp — i 



■—G R ^- I* 1 



6 

5a e 
5 



5() '' l " vA 
6 



((\a c ) = $dcvA(C|ae), 



6a' Gr JC )*dcvs(C'K) = $dcv£i(C'K), 



(64) 



cf. Eq. (j3"Tj) , and we arrive at Eq. (|ij9"|) with A c — > a c . 
This proves that Eq. (|53l) equally holds for "dressed" 
quasiprobability distributions. 



IV. DISCUSSION: CASCADED SYSTEMS AND 
SUDARSHAN'S OPTICAL EQUIVALENCE 
THEOREM 

A. Generalised optical equivalence theorem 

Analyses in Sec. IIIII conclude the principal part of our 
investigation. They show that the formal response pic- 
ture, postulated in Refs. by mere analogy with the 
harmonic oscillator, leads to a physical response formu- 
lation of electromagnetic interactions of distinguishable 
devices. 

Moreover, quantum electrodynamics in response repre- 
sentation looks wholly classical. The only difference from 
classical statitistical electrodynamics is that the "quan- 
tum probabilities" p l and p are not bound to be nonneg- 
ative. Any relation for P-functionals in quantum electro- 
dynamics holds as a relation for probability distributions 
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FIG. 2: Typical cascaded arrangements: (a) photodetection; (b) photodetection with coherent amplification. Shading signifies 
self-action problems included into models of the devices. 



in classical statististical electrodynamics. This statement 
constitutes the optical equivalence theorem in its most 
general form. 

The inverse theorem should be formulated with cau- 
tion. Namely, any relation for probability distributions 
in classical statististical electrodynamics, which (i) does 
not rely on their nonnegativity and (ii) is consistent 
with quantum electrodynamics, holds as a relation for P- 
functionals in quantum electrodynamics. For an illustra- 
tion of this statement see Sec. IIV Fl below. 



accounted for exactly. 

For general cascaded systems, the chain of devices and 
fields may contain arbitrary number of links. In Fig. [2p, 
we depict schematically a typical cascaded system where 
a field emitted by a source interacts with a coherent quan- 
tum amplifier and is then detected. Quantum theory of 
the systems in Fig. [2] is the subject of this section. 



Generalised photodetection theory without the 
rotating wave approximation 



B. Cascaded systems 



The model 



Taken that far, the optical equivalence theorem is of 
fundamental importance but of little practicality. For an 
arbitrary pair of devices, solving Eq. (|59[) is hardly easier 
than solving the dressing relation (f5T|) for the composite 
device. Things change if the macroscopic electromagnetic 
interaction becomes directional. Consider, for instance, 
the typical photodetection arrangement (Fig. [5^). Com- 
pared to the model of Sec. IIIIl this arrangement is subject 
to two additional approximations. Firstly, that one may 
distinguish the high-frequency (HF) optical field from the 
low-frequency (LF) photovoltage, and, secondly, that ra- 
diation of the source does not depend on the presence of 
the detector. Strictly speaking, these conditions intro- 
duce the concepts of source and detector. 

Note that these conditions concern not only physics 
but also engineering. Physics implies that the situation 
is macroscopic (mesoscopic) : devices are distinguishable 
and separated by macroscopic distances, and all light 
beams may be controlled. Proper engineering takes care 
of such problems as the influence of light reflected off the 
detector input window on the laser source. 

The photodetection arrangement in Fig. [2^ is a partic- 
ular example of a cascaded quantum system [12j ]. In pho- 
todetection, there are two distinguishable devices, source 
and detector, interacting with two distinguishable elec- 
tromagnetic fields, the input and output ones. Inter- 
actions of the devices with the fields are organised in a 
chain, which is made directional by neglecting the macro- 
scopic back-action of the detector on the source. All mi- 
croscopic electromagnetic self-actions are presumed to be 



What makes cascaded systems simple is that radiation 
incident on each link may be regarded given. Nontriv- 
ial self-action (dressing) problems are hidden inside the 
links, and may be approached separately. Assuming that 
the self-action problems for devices are solved, solution 
to a cascaded system reduces to a large extent to its ac- 
curate formulation. 

Consider, to start with, the generalised photodetection 
problem (Fig. [2ji) without the RWA. For simplicity, we 
treat the HF (input, i) and LF (output, o) fields as single 
modes, 



A(t) = 
A (t) 



2nh 
luoVo 
2irh 



(65) 



are 



H.c. 



where wo,i are the mode frequencies and Vo,i are the 
mode volumes. The field Hamiltonian reads, 



Hf — h(uioaQao +<Ji<zJai). 



(66) 



Generalisation to a more realistic case is straightforward. 

The model of Sec. IIIIl applies with minor amendments, 
due to the presence of two field operators. Hamiltoni- 
ans HdcvA.B an d the p-matrices PdcvA.B have the same 
meaning as in Sec. IIIIl One may regard i,o as two values 
of the variable x, so that 



J A (x, t) -> J A i(t), JAo(t), Jb{x, t) -t J Si (t), J Bo {t)- 



(67) 



9 



Device A interacts with the input field by means of the 
current J A i(t), while Jyi (i) — 0. Device B interacts with 
the input field by means of the current J B i(t), and with 
the output field by means of J B o(t) = Jo(t)- 

To unify the bookkeeping we introduce a Hamiltonian 
with "jumpers," 

H(t) =H { + H dcv (t) - [Ai(t) + Mt)] Ji(t) 

- [A {t) + A (t)]J (t), (68) 

where 

Hdev(t) = s A H dcvA (t) + s B H devB (t), 

Ji(t) = s A J Ai (t) + s B j Bi (t). (69) 

Jo(t) = S B J Bo (t). 

The "jumpers" sa.b =0,1 "commute" the problems. For 
example, with sa — 1, s B = we recover the problem of 
a solitary source. In fact we have to distinguish the prob- 
lem of device A and that of the light source, which differ 
in whether the c- number source Ai(t) is nonzero or zero. 
Correspondingly we have to define two types of quantities 
(averages): with nonzero A[(t), denoted (• • • )a, and with 
zero Ai(t), denoted (• • • )„. Specifications at the averages 
apply in fact to averaged operators, while quantum av- 
eraging as such is always over the p-matrix, 

p = |0)(0| g) /5 dcv = |0)(0| ® p AcvA ® PdcvB- (70) 

Redundant degrees of freedom are traced out automati- 
cally. For a summary of all definitions see table HI 

Full formal analysis of this system is the subject of 
appendix [Cl For purposes of this discussion, the source is 
described by the time-normal averages of the Heisenberg 
field operator Ai(t), 

(T :Ai(ti)---A i (t m y.} s 

=n{/ ^i(«)}p s (A)A(«i)---^i(^), (71) 

where we have introduced the corresponding quasiprob- 
ability distribution. The detector and the full arrange- 
ment are both described by the time-normal averages of 
the Heisenberg current operator j {t), but defined under 
different conditions. For the detector, 

(T :Jo(ti)---Jo(t m ):) d 

= II { / <^o(t)jpi(J \Ai)j (ti) ■ ■ ■ J {t m ), (72) 

while for the photodetection arrangement, 

(T ■Jo(tl)---Jo(trn):) 

= II { / d Jo{t)}po(Jo)Mh) ■ ■ ■ Mt m ). (73) 



The symbols {■ • • )d an d (■ • • ) D are defined in table|TJ Av- 
erages (|72)l by construction depend (are conditional) on 
the source Ai(t); this dependence is made explicit in the 
conditional quasiprobability distribution Pd{J \Ai) char- 
acterising the detector. Averages (|7T1) and (|73"1) are un- 
conditional. 



2. Semiclassical versus quantum photodetection theory 

Following the pattern of "doing quantum electrody- 
namics while thinking classically," the relation connect- 
ing observable photodetection statistics to properties of 
the source and detector should be, 

Po(Jo) = II { / ^i(*)jft(^i)Pd(Jo|A)- (74) 

Taking notice of the definitions ([71]). ([72]) and ([75]). Eq. 
(1741 is equivalent to the following relation for the time- 
normal operator averages, 

(T:J (h) ■ ■ ■ Jo(t m )i) 

= (T: [(T-Joih) ■ ■ ■ Jo(tm):)\ d ] \ A ^ A A- (75) 

For justification of these formulae see appendix [Cl 

It is instructive to compare Eq. ([75)) to the correspond- 
ing classical formula. In classical stochastic electrody- 
namics, a photodetector may be characterised by aver- 
ages of the photocurrent J(t) conditional on the detected 
field Ai(t), 

J(ti)---J(t m )\ [Ai y (76) 

In turn, the detected field is characterised by the uncon- 
ditional averages, 

Ai(ti)---Ai(t m ), (77) 

The unconditional photodetection statistics observed in 
the experiment is found imposing the second layer of av- 
eraging over the conditional averages ([76[) 

A. 

J(tl) ■ ■ ■ J(t m ) = J(tl)---J(trn)\ [Ai] '■ (78) 

Variable at the end of the bar is the one over which the 
averaging is performed; such specifications help to avoid 
confusion. Parallelism between Eqs. ([75]) and ([78]) is in- 
disputable: they coincide up to the replacement of time- 
normal averages by classical statistical averages. 

The said parallelism gets even more pronounced if we 
note that photocurrent must be classical also in quan- 
tum theory. That is, there must exist representations of 
the time-normal current averages by classical statistical 
averages, 

(T :J (ti)Jo(t m ):) d = J{h) ■ ■ ■ J{t m )\ [Ai]J (79) 
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Problem 


"Jumper" settings 


Relevant 


Notation for 


averages 






ext. souces 


"Raw" 


"Physical" 


Light source 


SA = I, SB = 


(Ai) 


(■■■U 


(••■ 


) 8 = ((-">^)| Ai=0 


Detector 


SA = 0, S B = 1 


A u (A Q ) 


(-)b 


(••■ 


)d = ({■■■)b)\ Ao=0 


Photodet. arrangement 


SA = 1, Sfl = 1 


(A,A ) 


(...) 


(... 


)o = (■ ■ ■ )| A . =Ao=0 



TABLE I: Three problems relevant to the arrangement in Fig. "Raw" averages imply the density matrix (|70|l and Hamil- 
tonian (|68|) . the latter with "jumpers" set to listed values. "Physical" averages follow by setting some or all c-number sources 
to zero. The table also lists the c-number sources on which the "raw" averages depend; those shown in brackets are set to zero 
in "physical" averages. 



and 

(T:Jo(tx)Jo(tm)-) = J(h)---J(t m ). (80) 

The same may be expressed as nonnegativity of the cor- 
responding P-functionals. Classical interpretation of the 
quantum averages ([79"]) and ([SU]) is possible if 

Po(J)>0, (81) 

and 

Pd(j\A)>0. (82) 

Conditions (IBTj) and ([52"]) are a heuristic principle to be 
imposed on quantum models of photodetectors. Nothing 
in the above and in appendix [C] depends on them. 

The only source of quantum behaviour in a photode- 
tection experiment is possible nonpositivity of the P- 
functional related to averages of the detected field, 

(r :M(ti)---Ai(t m y.) s 

= n{/ ' dA i (t)j Pa (A i )A i (t 1 )---A i (t m ). (83) 

Accounting for Eqs. ([79]), ([80]), Eq. (75]) may be written 

as, 

J(ti) ■ • • J(t m ) = | J dA { {t)^ Pa (A^ 

x J(h)---J(t m )\ [Ai] . (84) 

If p s (Aj) > 0, functional integration over p s [Ai) may be 
interpreted as a classical statistical averaging, 

(T:A{h) ■ ■ ■ A{t m ):) s = A i {t 1 )---A l {t m ) Al . (85) 

In this case Eq. ([84]) reverts to the classical formula ([78]) , 
and quantum mechanics becomes fully hidden from view. 

Noteworthy is that condition ([8"2"]) does not warrant 
condition (|8ip . Photodctcction must add enough noise 
to counter possible nonpositiveness of p s (A{). This is 
the reason why a photodetector free of shot noise cannot 
exist. 



D. Quantum and classical devices 

It is convenient to introduce terminology generalising 
the concepts of classical and quantum states of the har- 
monic oscillator to arbitrary electromagnetic devices. A 
device is macro scopically classical if the conditional P- 
functional characterising it is nonnegative, otherwise it 
is macroscopically quantum. Unlike the harmonic oscil- 
lator, macroscopic classicality or quantumness of a device 
concerns not only its quantum state but also its dynam- 
ics. We have to know both the state of the device pdev 
and the bare current operator J(x, t) in order to decide 
to which variety the device belongs. This knowledge suf- 
fices: the dressing formula ([25]) and Eqs. (O, ([52]) de- 
scribing interactions of devices preserve nonnegativity of 
P-functionals (i.e., if p l > 0, then p > 0; if p\ B > 0, 
then p l > and p > 0). 

The concept of macroscopically classical device calls 
for a word of caution. Firstly, any such device is quan- 
tum at a deeper level of insight. Classicality is always an 
approximation. Secondly, it may take quantum mechan- 
ics to explain properties of a macroscopically classical 
device. The best known example is the black-body radi- 
ation: by itself, it is in a classical state. Thirdly, there ex- 
ist devices in classical statitistical electrodynamics which 
cannot be implemented as macroscopically classical de- 
vices in quantum electrodynamics. Again, examples of 
such devices are well known: a coherent quantum ampli- 
fier, a photodetector without shot noise, and a light-beam 
doner, to name just a few. 

One more reservation one has to make is that a device 
may appear macroscopically classical due to limitations 
of an experiment. Indeed, let us have a closer look at 
functional pd(^o|^i)- It is natural to say that the de- 
tected field is in a coherent state if 

Ps (Ai) = l[6{Mt)-Mt)), (86) 
t 

where A e (t) is a given c-number field. For such field, 

Po(Jo) =p d (J \A c ). (87) 

That is, Pd(Jo\A c ) describes results of a photodetection 
experiment with light in a coherent state 18]. The pos- 
itivity condition (|82p warrants that the photocurrent in 
such experiment stays classical. However, Pd(^o|^4o) does 
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not give a full quantum description of the detector as a 
quantum device. Positivity of p&{J\E, E*) may be a re- 
sult of forfeiting all information about the optical mode, 
and/or ignoring response properties of the photocurrent 
mode (cf . appendix IC 2[) . The same device may behave 
classically in one experiment and in a quantum manner 
in another, depending on which information a particular 
experiment allows access to. 

E. Sudarshan's optical equivalence theorem 

That Eq. (|84p equally holds in classical mechanics with 
nonnegative p s (Aj) and in quantum mechanics with (pos- 
sibly) nonpositive p s (Aj) is Sudarshan's optical equiva- 
lence theorem generalised to interacting electromagnetic 
field. In Sudarshan's seminal paper [l9j], the optical 
equivalence theorem is formulated as follows. Let X(-, ■) 
be a normal representation of an operator X, 

X = :X{a,a i ):, (88) 

and P(-, •) — the diagonal representation (P-function) of 
the rho-matrix p, 

p = J d 2 zP{z,z*)\z){zl (89) 

where \z) is a coherent state, 

\z) = cxp(za t -z*a)\0), (90) 
and |0) is the vacuum state, 

a\0) = 0. (91) 

The quantum average of the operator X may then be 
written as a classically looking phase-space average, 

(X) = TrpX = J d 2 zP(z lZ *)X(z,z*). (92) 

Equation (|92p is a somewhat modernised version of eq. 
(5) in Ref. [l{|, which was written for a special case X = 
(v x a^. Parallelism between Eqs. ([51)) and (1U21 is evident. 

Sudarshan's optical equivalence theorem for interact- 
ing fields equally applies to classical and quantum de- 
vices. Indeed, Eq. ([71)) does not rely on heusistic con- 
ditions (f8"Tj) . (|52"|) . Equation ([75)1 may be written in yet 
another form, 

(r :jo{ti)---jo(t m y.) o 

= Il{/ d Mt)^Ps(Ai)(T:J (ti) ■ ■ ■ Jo{t m ):)\ d - (93) 

This relation holds irrespective of whether the time- 
normal averages here may be interpreted classically ac- 
cording to Eqs. (|79l) . (1801) . Any macroscopic device, no 
matter whether quantum or classical, is fully specified by 
its response to input field in classical states. 



F. Semiclassical theory of photodetection with 
coherent quantum amplification 

1. Preliminary remarks 

Utility of classical devices is in that they afford a phe- 
nomenological classical description within quantum elec- 
trodynamics. For a photodetector, this description is 
given by conditional statistical averages (I79|) . or, which 
is the same, by the functional p<\(J \Ai). It is often possi- 
ble to construct a simple phenomenological model which 
correctly reproduces averages (|79p . or analogous quanti- 
ties in other situations. If two or more classical devices 
are combined in a set-up with some quantum devices, 
a classical approach to the classical part of the set-up 
may preceed a fully quantum treatment of the set-up as 
a whole without any loss in rigour. 

As an example we consider the cascaded system in Fig. 
[2jx Subject to self-noise of the amplifier being in a clas- 
sical state, both the amplifier and detector are classical 
devices: if their input signal is in a classical state, their 
output signal is also in a classical state. Since, according 
to the optical equivalence theorem, full quantum descrip- 
tion of a device is given by its reaction to classical fields, 
theory of the composite detector (cf. Fig. [5b) may be for- 
mulated in purely classical terms, which is next to trivial. 

We also take this opportunity to generalise our analy- 
ses to the rotating wave approximation for optical fields. 
Consequently the HF fields and 1 (cf. Fig. [5b) are de- 
scribed by slow amplitudes E$(t) and Ei(t), which in 
quantum treatment become Heisenberg operators £o(t) 
and £\{t). The photocurrent J(t), which in quantum 
treatment is represented by the operator j7o(i), is a 
broad-band process and thus in no way is a subject to 
rotating wave approximation. 



2. Semiclassical versus quantum theory 

Full quantum treatment of the system in Fig. ^jp may 
be found in appendix [D] The result is very much as ex- 
pected: "doing quantum electrodynamics while thinking 
classically" is equally applicable to this system. In classi- 
cal stochastic electrodynamics, the detector is described 
by the probability distribution of the photocurrent, con- 
ditional on the detected field Ei(t), pd(J \Ei,Ei). The 
amplifier is described by the probability distribution of 
the output field, conditional on the input field Eo(t), 
p&(Ei, E* \E , Eq). As to the source, it is described by 
the unconditional probability distribution of the radiated 
field, p s (E , Eq) . It is convenient to start from solving for 
the composite detector (cf. Fig.[5j3), characterised by the 
probability distribution of the photocurrent, conditional 
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on the detected field E (t), p c d(</o|-Eo, Eq). Obviously, 



3. Classical coherent amplifier with spontaneous noise 



p cd (Jo\E ,E*) =Il{/ 

xp a (E u El | E , E*)p d ( J G | ^ , ) . (94) 

Statistical properties of the photocurrent, characterised 
by the unconditional probability distribution p (j ), 
then follow with ease, 

Po(Jo) = Il{/ d 2 E (t)^p s (E ,E*) Pcd (j o \E ,E*). 

(95) 

Alternatively, one may solve for the composite source (cf. 
Fig. Hb), characterised by the unconditional probability 
distribution p cs (Ei ,Ef\, 



Pc.(Ei,Et)=Y[{J d 2 E 



(t) 



and, 



Po(Jo) 



n 



x Ps (Eo , E*) Pa (E 1 , E\ \Ev,E%), (96) 



d 2 Ei(t) 



Pcsfa, E*)p d (j \E!, EZ) . (97) 



Either way, we find, 



Po(Jo) = H[j d 2 E (t)d 2 E 1 (t) 

x Ps (E , E*)p a (Ei , El | E Q , E*) Pd (J \E U E{ ] 



(98) 



In appendix [D] we show that Eqs. (J94|) - (j98|) equally 
hold as relations for P-functionals. A major fraction of 
the effort is spent on proper quantum definitions. How- 
ever, as soon as we become aware of this result, 90% 
of the quantum treatment turns redundant. Subject to 
self-noise of the amplifier being in a classical state, both 
amplifier and detector are classical devices. Their prop- 
erties may be described within simple semiclassical mod- 
els. The only possible source of quantum behaviour is the 
light source. To include it in our analyses, it suffices to 
assume that there exists a quantum framework, in which 
one can define time-normal averages of £o(t), 



<r:£b(ti)--.fib(*m)£o t (*i)---£o t (00 



(99) 



Under the rotating wave approximation, one has to apply 
apply the Kelley-Kleiner definition [t| [l(| , cf. appendix 
IA 31 Equation (|99[) is the only quantum formula we need 
in this section. 



The system in Fig. [2}d is rich enough in the sense that 
it allows one to illustrate restrictions to "doing quan- 
tum electrodynamics while thinking classically." These 
restrictions are two: quantum mechanics does not allow 
for noiseless amplification nor for noiseless detection. So, 
the simplest semiclassical model of the coherent amplifier 
consistent with quantum mechanics reads, 



Ei(t) = y/f~ a E (t) + E a (t), 



(100) 



where T a is the transfer coefficient and E a (t) is the ran- 
dom radiation (noise) added by the amplifier. This noise 
does not depend on Eo(t). If T a < 1 (attenuator), the 
noise may be neglected without getting into contradiction 
with quantum mechanics. (This does not mean that for 
T a < 1 the noise is always negligible: e.g., active medium 
without inversion attenuates the signal and adds spon- 
taneous noise.) If T a > 1 (amplifier proper), there is a 
well known limit pol ] to how small the noise may be (dis- 
cussed in Sec. IIV F~8l below) . This limit cannot be found 
classically, but there is no problem with including it into 
the classical Eq. (jlOOl) . This is the reason we call the 
model (|100p not classical but semiclassical. It operates 
with classical notions while making provisions for later 
"promotion" to quantum mechanics. By itself, Eq. (|100[) 
is fully classical, and may be handled without any refer- 
ence to quantum mechanics. 

For E a (t), we assume Gaussian statistics specified by 
the averages, 



E a (t) = 0, E a (t)E a (t') = 0, 
E* a {t)E a (t') = 27r/ a e-^/ 2 >l t - t 'l 



(101) 



We treat / a , j a and T a as phenomenological constants, 
available, e.g., from experiments with the amplifier. For 
consistency we have to assume that j a is large com- 
pared to the characteristic spectral width of the ampli- 
fied signal (otherwise one cannot regard T a as frequency- 
independent). That is, for our purposes, 

(102) 



E*(t)E a (t')=27rI aw 6(t-t'). 



where 



44 

7a 



(103) 



is the spectral density of the noise. Of use will also be 
the average, (recall that the noise is assumed Gaussian) 



\E a (t)\z\E a (t>W = \E a (t)\* 
= 4 7 r 2 / a 2 (l + e-^l*- t 'l) ~ 



\E a (t>)\ 2 
4^ 2 / a 2 [l- 



E*(t)E a (t') 



27a 



^(t-t')]. 
(104) 



4- Classical photodetector with shot noise 

In the conventional semiclassical theory of ideal pho- 
todetection [ll| , the inherently quantum nature of detec- 
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tion process is accounted for phenomenologically, postu- 
lating that photocurrent is a sequence of infinitesimally 
short photodetection pulses. The probability for a pho- 
todctection pulse to arrive within a sufficiently small time 
interval t, t + At equals 



p(t,t + At) = x \E 1 (t)\ 2 At. 



Here, 



X 



(105) 



(106) 



where rj is the detector efficiency and S is the beam area 
(assumed smaller than the detection surface). Planck's 
constant in (|106l) is all that survives of the quantum na- 
ture of the detector. If pit, t + At) is not small but E\ (t) 
may still be regarded constant on the scale of At, pho- 
todetection pulses within At obey Poissonian statistics. 
The bottom line is, all photodetection events (pulses) 
are independent of each other, and all correlations in the 
pulse sequence may only be due to correlations in the 
detected beam. 



5. Photocurrent statistics 

Based on these assumptions we can calculate all sta- 
tistical averages of the photocurrent. For averages con- 
ditional on Ei(t) we have, 



J (%E u E 1] =X\El(t) 



q- 2 J{t)J{t')\ [EuEl] = X \E 1 (t)\ 2 S(t-t') 

+ X 2 \E 1 (t)\ 2 \E 1 (t / )\ 2 , 



(107) 



etc., where q is the charge in a pulse. Using Eq. (|100|) and 
averaging over Eq{€) and E a (t) we find the unconditional 
(observed) photocurrent averages, 



J(t) =r] I T a n + n a 



2 J(t)J(t')=r } (T a n +n a )S(t-t') 

+ V 2 {T 2 n 2 (t - t') + 2T a n 0? i a [l + 4 7 " 1 S(t - t')] 

+ nl[l + 2 7 - 1 S(t-t / )}}, 



(108) 



where we have introduced the notation, 



S\E (t)\i 
n Q = , n 



S\Eg(t)f 

2irhuJo 



n 2 (t-t') 



S 2 \E a (t)\z\E (t')\2 



(109) 



For simplicity we assumed that E$(t) is a stationary ran- 
dom process. 



From a quantum perspective, hq and n a are the initial 
and added-noise photon fluxes and n^t — t') are autocor- 
relations in the former. Strictly speaking, at this stage 
we are not allowed to use this terminology, but we for- 
feit full rigour in favour of clarity. Another incidence of 
"premature quantum awareness" is separation of rj from 
h in Eqs. (|108|) . (|109[) . The parameter characterising the 
detector semiclassically is their ratio h/r]. 



6. Full classical solution 

General formulae for photocurrent averages are rather 
tangled, but a solution in terms of characteristic func- 
tionals is constructed with ease. Using the well-known 
characteristic function of the Poissonian process [2l| , for 
the functional <&&((,\E±,E{), replacing (|C14[) under the 
RWA for the detected field, we obtain, 

* d (C|£i.£i) =exp J dix|£iW| 2 { ex PfeCW]-l}, 

(110) 

It is equally straightforward to obtain photocurrent av- 
erages conditional on Eo(t). Formally, this corresponds 
to regarding the amplifier and detector as a composite 
detector (cf. Fig. [2Jd). A generating functional of such 
averages follows by making substution (|100[) in (|1 10[) and 
averaging over E a (t), namely, 



$ cd ((\E Q ,E*) 



= cxp J dtx\VT~*E (t) + E a (t)\ 2 {ex V [q((t)} -1} . 

(Ill) 

Observed photocurent statistics follows by averaging over 
the statistics of radiation of the source, 



3o(C) = $cd(C|£o,£ *) 



E 



(112) 



7. "Quantisation" of the problem 



"Upgrade" of the semiclassical theory to quantum 
fields is trivial. The composite detector (cf. Fig. [5}d) is 
a classical device, and the functional (1111[) may be "im- 
ported" to quantum electrodynamics as is. Proof of this 
statement, which takes quantum formulation of the ar- 
rangement in Fig. O), may be found in appendix [Dl (see 
also Sec. IIVF 2[) . All we need is to redefine averages 
(|112l) . namely, 



$ G (C) =(r:$cd(c|^o,4 t )=) s 



(113) 



The Heisenberg field operator So (t) is defined for the soli- 
tary source interacting with a narrow-band field. Hence, 
for purposes of this discussion, it suffices to assume that 
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there exists a quantum formulation of the source, assign- 
ing meaning to averages (j9"9l and thus to Eq. (|113[) . For 
formal particulars see appendix [Pi 

For Eqs. (|108j) . "quantisation of the source" reduces to 
redefinition of quantititcs (| 109|) : 



n 



S(£ f Q (t)£ (t)) s 



nl{t-t') 



S\TdlWl{t'%(t>%{t):) a 
^ 2 h 2 uj 2 



(114) 



We omit the symbol of time-normal ordering where it is 
redundant. If necessary, averages of E a (t) may also be 
redefined, e.g., 



S[(£l(t)Si(t)) 



aJ |£ =0 



(115) 



The average here is for a solitary amplifier without the in- 
put signal; as a formal quantity it is defined in appendix 
ID] While aesthetically satisfying, such redefinition is ut- 
terly unnecessary (so far as the self-noise of the amplifier 
is in a classical state). 

Alternatively, we can find photocurrent as dependent 
on radiation of the composite source (cf. Fig. [2b) , 

*«(C) 



Trexp / dtx£l(t)£i(t){exp[q((t)]-l}: 



(116) 



The symbol (• • ■ ) cs is defined regarding the source and 
detector and the electromagnetic interaction between 
them as a composite source (cf. Fig. [2Jd). For formal 
particulars see appendix [D] From Eq. (I116[) relation we 
find photocurrent statistics in terms of radiation of the 
composite source, 

q - 1 7ttj = x (£i(t)£ 1 (t)) ca , 

q - 2 J(t)J(t>) = x {£l(t)£i(t))J(t - f ) 

+ x\T^£\{t)£\{t')£ l {t l )£ l {ty.) ca . 

(117) 



As expected, Eqs. (|117[) are standard Glauber-Kelley- 
Kleiner's photodetection formulae [ill ]. 



8. Limits on spontaneous noise 

Quantum mechanics imposes a lower limit on the spec- 
tral density of spontaneous noise, 



4n a 1 r - \ 

/a ^ 



T n -1 




For amplifiers based on stimulated emission, this is a con- 
sequence of the connection between the Einstein coeff- 
cients, cf. eq. (44.9) in Landau and Lifshitz's textbook 
[22| . As a general quantum relation, this formula was 
deriv ed, e.g., by Caves (specifically, this is eq. (4.21) of 
Ref. |2(| rewritten in our notation, cf. endnote |23|). 

To see how condition (|118[) serves to keep the semiclas- 
sical theory consistent, consider the spectrum of fluctua- 
tions of the photocurrent, 



J 2 = / rfre MT J{t)J(t + r). 



(119) 



Using Eq. (|108l) and ignoring the zero-frequency (con- 
stant) contribution we obtain, 

q~ 2j u = v( T ^ n o + n a) 

+ V 2 {T 2 n 2 0LU + 8T a n n a7a - 1 + 2n a7a x ), w > 0, (120) 
where 

S ' r dr ^ (T:£l ml (t')£ (t')£ (*):>. 



(121) 



Unlike in classical mechanics, in quantum mechanics Uq^ 
is not bound to be positive. There is an obvious limit to 
how negative this quantity may be. Without the ampli- 
fier, photocurrent spectrum reads, 



q 2 J2 = vno + rf n 



0u> ■ 



(122) 



Photocurrent is classical, so that its spectrum should be 
positive. Assuming an ideal photodetector, r\ — 1, we 



find the lower limit on riQ U , 



n 2 0uj > -n . 



(123) 



Furthermore, assume that at u> = u) this limit is actually 
reached, 



n oa> = -no- 

At this frequency, Eq. (|120[) yields, 



(124) 



Q 2j l = V T a. n o(^ - + 8^a7a ^ 

+ r;n a (l + 2r / n a7a - 1 ). (125) 

Having this quantity positive for arbitrary no and 
r] < 1 imposes a lower limit on the spectral density of 
spontaneous noise, 



4n a 1 / s 



(126) 



This condition is weaker than (|118l) , so that amplification 
is bound to degrade quantum properties of the signal. 
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V. CONCLUSION 

So what is it that we wish to say and that, we be- 
lieve, has not been fully appreciated for more than 40 
years? Conventional wisdom is that relation like Eq. 
(jllOP are limited to classical input fields. This is in- 
deed the case if we regard it as characterisation of the 
photocurrent . However, regarded as characterisation of 
the photodetector , Eq. (|110[) is a wholesome quantum for- 
mula. Moreover, it is a complete quantum formula: ow- 
ing to Sudarshan's optical equivalence theorem, there is 
no freedom in generalising semiclassical photodctcction 
theory to quantum fields. It is this absence of freedom 
that we wish to stress. It applies to any device: if we 
know how it reacts to classical fields, we know how it 
will react to quantum fields. In principle, it applies also 
to quantum devices, but its most dramatic consequences 
are for classical devices. Sudarshan's optical equivalence 
theorem turns their semiclassical models into quantum 
approaches. 

In conclusion, quantum theory of the electromagnetic 
interaction under macroscopic conditions of distinguisha- 
bility of devices and of controlled actions and back- 
actions between them is constructed. This theory is 
subject to "doing quantum electrodynamics while think- 
ing classically," which allows one to substitute essentally 
classical considerations for quantum ones without any 
loss in generality. 



Appendix A: Summary of formal definitions 

1. Time orderings of operators 

By definition, bosonic operators commute under all 
kinds of ordering. The standard time-ordering, denoted 
T + , puts operators in the order of decreasing time argu- 
ments, 

T+X 1 (t 1 )X 2 (t 2 ) = e(t 1 ~t 2 )X 1 {t 1 )X 2 (t 2 ) + {l^2}, 

(Al) 

and similarly for more operators. The "reverse" time 
ordering, denoted T_ , may be introduced by the relation, 

[T+X^X^h) ■ ■ ■ X m (t m )f 

= T_xl{t 1 )Xl{t 2 )---Xl{t m ). (A2) 



The closed-time-loop, or C-contour, ordering [15Ml7l | is 
a way of writing the double time ordered operator struc- 
tures, 

Tcx^ih) ■ ■ • -Y m _(t m )j>i + (ti) • • • y n +(t' n ) 

= T_X 1 {t 1 )---X m {t m ) T+A {t[ )---X(t' n ). (A3) 



The ± indices in (|A3|) serve only for ordering purposes 
and otherwise should be disregarded. We note in passing 



that definitions like (|A1|) . (IA3|) may cause mathematical 
problems, see the concluding remark in appendix Al in 



2. The frequency-positive and negative parts 

The symbols ^ ' denote separation of the frequency- 
positive and negative parts of functions, 

f(t) = f i+) (t) + f(-\t), 

+oo +00 

/(±)(<) = / ^ e ~^( ±w )/- U = J die™* f{t). 

— oo — oo 

(A4) 

This operation is alternatively expressed as an integral 
transformation, 

/ (±) (*) = J dt'5^(t-t')f(t'), (A5) 

where 

^(t) = ^)(_ t)=[ ^) (t)] - = ± _l_^ (A6) 

are the frequency-positive and negative parts of the delta- 
function. For more details on this operation see Ref. [H, 
appendix A. 



3. Time- normal ordering 

Definition of the time-normal operator ordering, de- 
noted T: ■ ■ ■ :, is different for slow amplitudes, such as 
£{x,t) and V(x,t), and for Hermitian broad-band fields, 
such as A(x,t) and J{x,t) (for purposes of this discus- 
sion, their physical nature is irrelevant). To the for- 
mer, one applies the conventional definition of Kelley and 
Kleiner [H, [lC [_, while for the latter, the amended defini- 
tion of Refs. §, [|| must be used. Most cases of interest 
in the paper are covered postulating the operator-valued 
characteristic functional, 

T: exp (irjA + iQJ + i[i '£ — ifj£' + iv*T> — ivTr) : 
= T c exp [ir ] ( -- ) A + + + J+ + i( i+) J- 

+ i/i*£+ — in£_ + iv*T)+ — ivT>\ 
= T_ exp [i-q^A + iQ {+) J - iy£^ - ivV^] 

x T+ exp [irf^A + iC^J + i[i*£ + iv*f>\ , (A7) 

where r)(x,t), £(x,t), (x(x,t), v(x,t) are auxiliary c- 
number functions. We use here notation ([7]); separation 
of the frequency-positive and negative parts was defined 
in Sec. IA 21 Similar definitions apply to free operators; 
for E(x,t) and A(x,t) defined by (|B1[) the time-normal 
ordering reduces to the standard normal ordering. For 
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an in-depth discusion see Ref. Q. Causality properties 
of the time-normal ordering, which are of crucial impor- 
tance for physical consistency of our analyses, were the 
subject of Ref. [25|. 

Appendix B: The general case of interaction of two 
distinguishable devices 

1. The model 

In the general case, a subset of oscillator modes @ 
is made subject to the rotating wave approximation 
(RWA). All oscillators are organised in two quantised 
fields: the narrow-band, or resonant, field E(x,t), and 
the broad-band, or nonresonant, field A(x,t), (1 < M < 
AO 



&{x,t) = -i^^<{x)aiS^-^\ (Bl) 



A(x,t)= J^u K (x)a K e- iu ^ + H.c. . 

k=M+1 ' L ° K 

Frequencies ui K , 1 < K < M, are supposed to occupy a 
narrow band centered at uiq, so that E(x,t) is by defini- 
tion a slow amplitide. This assumption only matters for 
physics; formally, it may be disregarded. 

The general case is equally governed by the generic 
Hamiltonian (flj, where the electromagnetic interaction 
is now split into the narrow-band (RWA) and the broad- 



band (no-RWA) parts, 

Hi{t) = - j dx(\\A(x, t) + A c (x, t)] J[x, t) 

+ {[E(x,t) +E c {x,t)]D\x,t) +H.c.}). (B2) 

The Hamiltonian H dcv (t) in ([1]), the dipole momentum 
D(x,t) and the current operator J(x,t) describe the de- 
vice. They commute with all a K , aJ H and otherwise remain 
arbitrary. The c- number external sources E e (x,t) and 
A c (x,t) are added for formal purposes. For a discussion 
of this model see our Ref. @, sections II and III. Com- 
pared to Refs. [H, 0], interaction (|B2[) lacks the c- number 
dipole and current which are of no use in this paper. 

The broad-band field enters the theory through the 
retarded Green (response) function ([6]), while the narrow- 
band one — through another response function, 

A R {x,x',t-t') = ~6(t-t')[E(x,t),&(x',t')]. 

(B3) 

The operators E(x,t), A(x,t), D(x,t), and J(x,t) are 
by definition the interaction-picture (free) ones. Their 
Heisenberg counterparts will be dehoted as, respectively, 
£{x,t), A(x,t), V(x,t), and J{x,t). 

2. Quantum electrodynamics of a solitary device 
revisited 

Extension of formulae of Sec. |TT] to the general case 
reduces to a large extent to triplicating all variables. In 
place of definitions ([2~T|) . (|22|) we have, 



$*, ov (C, v, v*\A c ,E c , E*) = Trp dcv T: exp (itf' + W*D' - ivD*) : 

= n{/ dJ(x,t)d 2 D(x,t)^p I (j, D, D*\A e , E e , E*) exp (i(J + iu*D - ivD*), (B4) 



where the time- normal ordering is defined by Eq. (jA7[) . and the primed operators are defined as Heisenberg ones with 
respect to the Hamiltonian, 



H(t) = H dcv {t) - J dx[A c {x,t)J(x,t) + E*(x,t)D(x,t) + E c {x,t)& (x,t)]. (B5) 
This is Hamiltonian ([1]) with field operators set to zero. In place of (JTTj) we find, 
$dev((,v,v*\A e ,E e ,E*) = ^T:exp (i(J + iv*V - ivV*) :^ 

= II { / d J ( x ' t )d 2 D{x, t) \p(j, D, D*\A C ,E C , E*) exp (i(J + iv*D- ivD*) . (B6) 

x.t ^ •> ' 

The dressing formula reads 0, 

( S S 6 S S 5 \ 
$dev((,v,u*\a e ,e e ,el) = exp ^ - i 's a ~ GK S( ~ i Je~^ K ~5v* + i Je* ^Jv) $ dov(C' v , v*\a c , e c , e*), (B7) 
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with the equivalent formula for the quasiprobability distributions being Q, 

P (J,D,D*\A C ,E C ,E* C ) =p\j,D,D*\A c + G R J,E c + A R D,E: + A* R D*). (B8) 
Functionals $ x and <I>dcv may equally be defined with operators without sources. For the bare device, 

$dcv(C v -> v*\ a c e c , e*) = TrpdcvTc ex P («C+^+ - K-J- + iv + D + + iu + D\ - iv-D- - iv-D^_) I , (B9) 



cf. Eq. (|23p. where c.v. (short for causal variables) refers to the union of the nonresonant response substitution (|19j) 
and of the following resonant one, 

e c (x,t) _ e*(x,t) 

i/+(X,t) = - , ^+(2:, t) — V (X, tj H , 

• / S f\ ( B1 °) 



For the dressed device an extension of Eq. (|24|) may be derived, 

*dev(C> v -> v *\ a v + A c, e c + E c , e* + E*) 

= Tip dcv T c exp(i( + J + -iC-J-+iP + T) + +w + vl-iD^V^-i^vl)\ . (Bll) 

I C.V. 

For a verification of Eqs. (|B4|) - (|BTT|) see Refs. @, 3 0|- 

3. Interaction of distinguishable devices 

Conditions ([32 ]) -([3l ]) . ([35 ]) — ([37 )1 and (|38]) -(|50] l are supplemented by, respectively, D(t) = D A (x, t) + D B (x, t), 
D(t) = D A (x,t) and D{t) = D B {x,t). The bare and dressed functionals characterising the components of the device 
are defined inserting indices A, B into Eqs. ()B4| - (jB8|) . The connection between p 1 and p\ B is a natural generalisation 
of Eq. p]), 



o l (j,E,E*\A u D h D?) =Jj| / dJ A (x,t)dJ B (x,t)d 2 D A (x,t)d 2 D B {x,t) 



x <5(j(>,t) - J^OM) - J B (:M))(5 (2) (L>0M) - D A {x,t) -D B (x,t)) 

x p^(j A , D A) £> A | Ai, Ei,E^)^ B {J B ,D Bi D* B \A h E h Ef). (B12) 
Using the general dressing formula (|B8[) and proceeding as in Sec. IIIIBI we arrive at a generalisation of Eq. (|53p . 

p(j,D,D*|A e! D e! D*) = Jjj f d.J A (x,t)dJ B {x,t)d 2 D A {x,t)d 2 D B (x,t) 

x S(j(x,t) - J A (x,t) - J B (x,t))S {2) (D(x,t) - D A (x,t) -D B (x,t))^ 

xp A (j A ,D A , D* A \A cA , E cA , E* cA )p B (J b ,D b ,D* b \ A cB , E cB , E* B ) , (B13) 
where A cA ,b are given by Eqs. (|52|) . and E eA ^ B — by the analogous formulae, 

E cA (x,t) = E c (x,t) + / dx'dt'A K (x,x',t-t')D B (x,t'), 



(B14) 

E eB (x,t) =E e (x,t) + I dx'dt'A R (x,x',t-t')D A (x,t'). 



In operator terms Eqs. (|B12[) and (|B13p are equivalent to, correspondingly, 

$ dev (C> v, v*\a Cl e c , e*) = $ dcVj4 (C, f, ^* \a c , e e , e*) $^ ovB (C, f, v*\a e , e , e*) , (B15) 
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and 

^ /> *l ^ ( ■ 6 A S ■ S A 6 ■ S A* 5 ■ S A* S 

$ dcv (C, „, v \a e , e c , e c ) = cxp ^ - ^A R — - ^A R - + ^A R — + ^A R - 

- i- — G-R-fr: - iJ-rG n -^-] <S> dcvA ((i v, v* \a e , e e , e*)$ dc vi3 (C, v'* \a' e , e' e , e' c *) I , (B16) 

da e dQ' da' c Of, J 1 £'=<y=v,a'e=a<>,e-i=ee 

I 



cf. Eqs. (135)1 and (P)) . 

Appendix C: Photodetection problem without the 
rotating wave approximation 

1. Formal quantum solution 

In this appendix we give a formal justification to Eqs. 
(|74|) and ([75]) . In physical terms, this demonstrates con- 
sistency of "doing quantum electrodynamics while think- 
ing classically" with the approximation of avoided macro- 
scopic back-action of the detector on the source, which 
underlies the very concept of photodetection. For all def- 
initions see Sec. II V Cll and table |U 

We adapt the algebra of Sec. IIIIDI The presence of 
two broad-band fields is formally accounted for assign- 
ing indices i,o to the field and current operators and all 
auxiliary quantities, and redefining condensed notation 
accordingly, e.g., 

CJ= /dt[G(t)ii(t) + Co(*)£(t)], (Ci) 

etc. Such redefinitions automatically extend all results 
of Sec. IIIII for the broad-band case to the problem at 
hand. The emerging formal structure is anything but 
transparent, and, worse, takes no heed of specifics of the 
problem at hand. We therefore keep indices i,o explicit. 
Bare devices are characterised by the functionals, 

^(Ci|A) =(T:exp(iCiJ[):) A , 
$dcvij(Ci,Co|^i,A ) = (T:exp (iCiJ[+iCoJo)-) B , 
$ dcv (Ci,Co|Ai,A ) = (T:ex.p(i{ i j! + i(J' ):) 

= *dcvA (Ci | A) $dcvB (Ci , Co | A , Ao) , 

(C2) 

where the primed operators are the Heisenberg ones for 
the Hamiltonian dHHJ) with A^t) = A a (t) = 0. Device 
A does not interact with the output field; consequently 
3> deVj4 does not depend on £ (i) nor on A Q (t). 

For the dressed devices, use will be made of the func- 
tionals, 

$dcv4Ci|A) = (T:exp(iCiJi):) A , 
$devB((u(o\Ai,A ) = (T :exp (i^Ji + KoX)-) B , 
*dcv(Ci,Co|Ai,A ) = (T :exp (KiJi + KoX)-)- 

(C3) 



Calligraphic letters denote Heisenberg operators. These 
are in essence placeholders, the exact meaning of which 
depends on "jumper" settings (see Sec. IIVC"T|) . 

Response function Gr given by © is now a 2 x 2 matrix 
of kernels, however, since A{(t) and A Q (t) by construction 
commute, this matrix is diagonal. All quadratic forms 
involving Gr split in two, 

rjGnj = ViGmji + ?7oG Ro j , (C4) 

and the causal reordering exponent factorises, 




The kernels Gr; and Gr are defined applying (|6]) to the 
input and output fields, 

G m {t-t')= l -e{t-t')[A l {t),A l {t% 

n (C6) 
G Ba (t-t') = ~6(t-t')[A (t),A (t')]. 

In place of Eq. (|60[) we thus have, 

$dev(Ci,Co|^i,A ) = 

° XP ( l lA GRi Jc) 6XP ( - l A- GR °^k) 

X $ dcvA (Ci|A)$dcvs(Ci,Co|A, Ao). (C7) 

The "output" exponent here acts only on $ devB ; Eqs. 
(15T]) - (p3")) should therefore be applied to the "input" ex- 
ponent (with C ^ Ci an d o-e — ► A{). The dressing rela- 
tions now read, 

6XP ( ~ 'JI^J^)^^^ = $dev ^(Ci| A )' 
exp ( - A) exp ( - i^ORo^) 

X $dovs(Ci',Co|4, A Q ) = $dovB(Ci',Co|^, A ), 

(C8) 

while Eq. (|59p expressing interactions of dressed devices 
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becomes, 

$dcv(Ci,Co|A,^o) 

exp 



.66 .66 
% 6^ Gm 6Z~ 6A m 6Q 



x $ dev A((l\A[)$ devB ((i,(o\A h A )\ (t _ CuA ,_ Ai . (C9) 



Approximations and solution to the 
photodetection problem 



Equation (|C9[) is not yet a photodetection formula. 
Firstly, the source-detector interaction in it is bi- 
directional. Secondly, it contains a lot of irrelevant in- 
formation, in particular, full quantum (response) prop- 
erties of the input and output helds. So, it "knows" how 
an attempt to measure the input held would affect the 
detection, and how simultaneous measurements of the in- 
put held and output current would be correlated. Subject 
to valid quantum models of the devices, it also "knows" 
about all limitations imposed by quantum mechanics on 
such simultaneous measurements. 

By definition, of interest to us is the output signal of 
the detector under the condition of avoided macroscopic 
back-action of the detector on the source. The said signal 
is formally expressed by the joint time-normal averages 
of the field A (t) and current J (t), defined in the ab- 
sence of external sources. These averages are accessible 
through their characteristic functional, 

(T: exp (ir) A Q + iC, J ) -) 

= $dcv(0,Co + ??oG Ro |0,0) =$o(Co + r/ G Ro ), (CIO) 

cf. Eq. (HU). The symbol (• • • ) D was defined in Sec. lIVCll 
(see table 

Equation (fCTOf suppresses the aforementioned irrel- 
evant information. The no-macroscopic-back-action ap- 
proximation is imposed dropping the signals propagating 
from device B to device A. Formally, this means a re- 
placement in Eq. (IC9|) . 



(Cll) 



• 6 „ 5 . 6 6 
ex Pl -i—G Ri -- l —G Ri - 

Under this assumption Eq. (IC9[) yields, 



(3 3 
~ 1 5A Gr % 

X $devA(Ci|0)$dovs(0,Co|A,0)| Ci=0 A[=0 . (C12) 



We dropped the prime at ([(t) which after setting 
Ci{t) = became redundant. 

Confining out interest to the output signal reduces the 
information required about the devices. The source is 



described by the time-normal current averages generated 
by the functional, 

$dev4Ci|0) = (T:exp(iQj A ):) s = (C13) 

while the detector is characterised by the time-normal 
current averages conditional on the input held, 

$dcvi3(0,Co|Ai,0) 

= (T:exp(< J Bo ):) d = $ d (Co|A). (C14) 



The symbols (• • • ) s and (• • • )d were defined in Sec. lIV C~T1 
(see table HJ). With irrelevant information hidden from 
view we find the photodetection formula, 



$o(Co) = exp ( - 



6 r 6 
6Q 



5Ai 



x$ s (Ci)^d(Co|A)| fi=0 , Al=0 . (Ci5) 



It is worth empasising that, to obtain Eq. (|C15|) . the 
"one-way" assumption does not suffice. We had also to 
assume what we do not do certain things in the experi- 
ment: do not attempt to perform any additional measure- 
ment on the optical mode, and do not attempt to probe 
the detector by radiation in the photocurrent mode. 



3. Photodetection statistics as a quantum average 
over the detected field 



It is straightforward to rewrite Eq. (IC15I) in terms of 
the detected field rather than the source current (which 
is a customary viewpoint). Time- normal averages of the 
held radiated by the source are generated by the func- 
tional, 

T s ( m ) = (T: exp (irjiAi) :) g = $ s ( m G Ri ) , (C16) 

where use was again made of Eq. (|15|) . Applying the 
obvious formula, 



$ s (?7iG R i) = exp ( ?7iG RiT7 - )$ s (Ci)|, =0 , 



1 6CJ 



we can rewrire (|C15[) as, 



*o(Co) 



SAi 



MCo\Ai) 



(C17) 



(C18) 



Application of the differential operator here cannot be 
anything but a fancy way of expressing quantum averag- 
ing over the detected field. Indeed, by definition 



^) = ( T:eXP ( ii ^ i 



(C19) 



cf. Eq. (IC16|) . Under orderings, operators behave as c- 
numbers, and the exponent in (|C19I) may be interpreted 
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as a functional shift operator. Hence for any c-number 
functional T(A{), 



[(T:r(Ai + A):)J | Ai=0 - <r=^(A):) s . (C20) 



Using this formula we can rewrite Eq. (|C18[) explicitly as 
a quantum average over the detected held, 



<fo(Co) = <T:$d(Co|ii):) s 



(C21) 



Equation (|75j) readily follows from this formula, which in 
turn gives justification to Eq. (|74p . 



by means of the dipole operator L>Ao(t)- Device B 
interacts with mode by means of D B o(t) and with mode 

1 by means of D B x(t). Device c interacts with mode 1 
by means of Dci(t) and with mode 1 by means of the 
current Jc2(t). The state of all oscillators is vacuum, so 
that the state of the full system is, 



p = |0)(0| ® /SdcvA ® /OdovB ® pdevC- 



(D3) 



Formally, we have to consider six physical problems: 
three of solitary devices, one of the composite source, one 
of the composite detector, and one of the whole system 
(cf. Fig. [2b)- To unify the bookkeeping we postulate the 
Hamiltonian in the form, 



Appendix D: Quantum theory of a cascaded system 

1. The model 

The high-frequency (HF) optical fields and 1 in Fig. 
[2]b are treated under the rotating wave approximation, 
while the photovoltage — which is a low-frequency (LF) 
field — without the rotating wave approximation. The 
corresponding held operators are, 



Eo(t) = i\— - — a , Ei[t) =«*/— - — oi, 



Vn 



A (t) 



2nh 
W2V2 



(Dl) 



a 2 



e~ tU2t + H.c. . 



where 60,1,2 ar e the standard annihilation operators, 
^0,1.2 are the frequencies and Vb,i,2 are the so-called mode 
volumes. Recall that 2?o,i(i) are slow amplitudes; they 
lack time exponents. The Hamiltonian of the electromag- 
netic field is a sum of three oscillator Hamiltonians, 



K = 



(D2) 



The Hamiltonians of bare devices are Hd ev A,B,c(t) and 
their states are PdevA,B,c- Device A interacts with mode 



H(t) = H S + H dcv (t) - [A (t) + A (t)] J (t) 

- { [E (t) + E (f )] L>1 (t) + [Ex (t) + Ex (t)] b\ (t) +H.c.} , 

(D4) 

where 

Hdev(t) = SAH<x ev A{t) + Sfi-ffdevB (0 + ScHdevC\t), 

D (t) = s A D AQ (t) + s B D m {t), 
Dx(t) = s B D m {t) + s c D C x(t), 

Jo(t)=S C J C 2(t). 

(D5) 

The "jumpers" sa,b,c = 0, 1 serve to "commute" the 
problems. For example, with sa = ^,s B — sc — we 
recover the problem of a solitary source. More precisely 
speaking, we have to distinguish the problem of device A 
and that of the light source, which differ in whether the 
c-number source Eo(t) is nonzero or zero. Correspond- 
ingly we have to define two types of quantities (averages): 
with nonzero Eo(t), denoted (• • • )a, and with zero Ea(t), 
denoted (• • • ) s . For a summary of all definitions see table 

m 



"Dressed" devices are characterised by the averages, 







= (T: exp 


{wqV 


- w Q t)l) 








[uo,Uq,vi,u*\Eo,Eq,Ei,E*) 


= (T: exp 




- ivoDl 


+ iv\Vx 


-iuxT>\):) B 




^deyc(^l,^xXo\Ex,El,Ao) 


= (T:exp 


(iutVx 


- ivxf>\ 


+ <o Jo) 




^devBc( J/ ) ^0' 


v\, v*, Co \Eq, Eq , Ex , E{ , A ) 


= (T: exp 


+ 


- WoDq 

iv\T>x — 


wxV>\ + 


i(oJo)'-) BC , 


$dov(^0, ^0, 


^i,^xiCo\E ,Eq,Ex,EI, A ) 


= (T: exp 


+ 


- %vqV\ 
iv\t>x ~ 


ivx±>\ + 


«Co Jo)-)- 
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Problem 



"Jumper" settings 



Relevant 
ext. souces 



Notation for averages 
"Raw" "Physical" 



Light source 


SA = 


= 1, SB - 


= 0,sc ~ 


= U 


(-Bo) 


■ ■ ■ )a 


1 

(■■ 


' )s — 


= ((■ ■ 


■ ) a >\e =o 


Amplifier 


SA = 


= 0,s B = 


= Mc = 


= 


E ,(E 1 ) 


■ ■ ■ )b 


(•• 


/a ~~ 


= ({■■ 




Detector 


SA = 


- 0, SB ~ 


= 0,s c = 


= 1 


E U (A ) 


■ ■ ■ )c 


(•• 


' )d = 


= ((•• 


■) c )|a o = 


Composite source 


SA = 


~ 1) Sb - 


= l,sc = 


= 


(B ,£l) 


■ ■ ■ )ab 


<•• 


' ) cs ~ 


= ((• 


' ) AB ) \e =e 1 =o 


Composite detector 


SA = 


~ 0, SB - 


= Mc = 


= 1 


So, (Ei, A ) 


■ ■ ■ )bc 


(•• 


• )cd 


= «• 


) Bct )|b 1= o,a =o 


Full system 


SA = 


- 1) SB - 


= l,sc = 


= 1 


(E ,E x ,A o ) 




<•• 


' /° ~~ 


= <••• 


'\En=Ei =0,A o =0 



TABLE II: Six problems relevant to the arrangement in Fig.[5]3. "Raw" averages imply the density matrix (|D3|I and Hamiltonian 
(|D4[) . the latter with "jumpers" set to listed values. "Physical" averages follow by setting some or all c- number sources to 
zero. The table also lists the c-number sources on which the "raw" averages depend; those shown in brackets are set to zero in 
"physical" averages. 



Calligraphic letters are as always for Heisenberg operators. For the "bare" devices, 



*dc V B (vo , vj) , "i , v\ I E Q , El , E l , E{ 



(T: exp (wq-Do 
(T: exp (iv%D' 
^\ cvC {u u ylXo\E u El,A ) = <T:exp (ii/JDj 

(T: exp (iv*D' 



^devsc (^o, ^o i V U v * > Co \Eq, Eq , Ei , E^ , A a ) 
*dev & , "a > v x , v\ , Co | E , E* , Ei , E{ , A ) 



iis D'J)--) A , 

ivoD'J + iv{D' x - iuiD^)\) B , 



+ iv*D[-iv 1 D? +i( Z):) 



BC 1 



(T: exp (wq-Dq - iu D^ 



KoJo)-)- 



(D7) 



The primed operators are the Heisenberg ones with respect to the Hamiltonian, 



H(t) = H f + H dcv (t) - A (t)J (t) - [E (t)E>l{t) + E^Dlit) + H.c] . 



(D8) 



In (|D6|) and (|D7|) . specifications at the averages apply in fact to averaged operators, while quantum averaging as such 
is always over the p- matrix (|D3[) . Redundant degrees of freedom are traced out automatically. 
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Formal solution 



Adapting the general dressing formula (|B7[) to the six problems at hand we have, 

^devA^O.^O \ E 0,Eq) = expZ ^^-, ^^devA^Qi^ol-^O,^), 



< i ) dovs(^o, Vq, i>i, v*\Eq,Eq, Ei,Ef) 
= exp 

^dovc o ,Vo,vi, v\ , Co | Ei , E{ , A Q 
= exp 

^devBC^Oi ^0' V l, U l ' Co|S(), Eq,Ex,EI,A 

= exp 



z M s \+z (— — 



*dcvB ("a , v*o , v\ , "I | S , E* ,E x ,El), 

*dovc(^^r,Co|^i,-Bi*,A ), 



Z ° ( M 'sl Q ) +Zl { A >sl) +Z °{ SA Q ' 6C Q 

X ^dcvBC («D . v > "1 . "l . Co | ^0 , S* . ^1 . E l > 4>) 
^dev(l / 0, ^O) V l, V l ) Co | -So j Sq, Si, £*, A ) 



exp 



Z 



5E(j ' Ji'o 



6 S \ (5 5 



$ Lv ("a , ^o i "i ) W i Co | S , Sq , Ei , El , A ) , 



(D9) 



where 



Zi 



6 _S_ 

5Eq ' Suq 
5 _S_ 
SEi ' Svi 



-°\5A >5Co 



A 



R0 



SE 

SEi ov{ 

■ 5 r 6 
- 1 6A Gr °6(-o- 



+ i 



-A 



SE* ™5v 
■ 5 a* 5 

l 6E* /Xm Su~i 



(D10) 



The kernels are given by the formulae, 



A R0 (t - 1 




O[So(t),s o t(t0], 


A m (t-t 






G Ro (t - 1 




0[ioW,io(0]- 



(Dll) 



Functionals (|D6I) may be redefined in terms of "bare" dipole and current operators without sources in the manner 
of Eqs. (|B9[) . (|B10[) . This leads to the factorisation properties, 



*dcvsc iyo , vl , vi , v\ , Co | S , El , E x , Et , Ao) = <t>\ cvB (vq , v* , v x , v\ | S , S* , E x , ) 

X ^devC^l.^l.ColSx.S^Ao), 
*dcv (MD , Vl , I/l , I/f , Co I S , -B * ,EuE*>A>) = ®d C vA (vq , vl I So , S*) 

x ^LvBC (" o. fo > " l . v* , Co I S , Sq* , Ei , E{ , A Q ) 
= ^ (i* , K | So , E*)^ dcvB (uo , ^ , vx , | So , S* , £ x , £* ) $^ cvC , „* , Co | Ei , i?* , A G ) . 



(D12) 
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Combining these formulae with the dressing ones (|D9I) we find the relations among properties of the "dressed" devices, 



$devsc ("o, , v x , v{ , Co I E , Eq , E x , E{ , A ) = exp 



8 Ei ' 8v[ 



Zl 777' 777 + Zl 777' 777 



x $ d cvi3 fa , , v 'i - "i* I . B o . > E'i ) $dcvc fa , v\ , Co | Ei , E{ , A ) I ^ =!/i £ , =£i , 



$dev(l / 0,I / 0) I/ l> I/ *)Co|-E , 0,-Eo>^l>- E; i>-^o) = exp 



Z 



8 _8_ 

8E '8i/ 



8 _S_ 

JeI'sVq 



X $dcvA fa , "0 I #Ch ^0* ) $dcvBC fa , ^ , Vl , V* , Co \Eq, #0* , ^1 , ^i* , 4>) | ^, = „ Q £ , =£q 



exp 



(5 5 
SE^'H 



SE' ' <Ji/„ 



g _5_ 



JW^sVi 



x$ devA (v' ,v' *\E' ,E' *) 

x $dcvB(^o,^o^^^^1^o,^o^^>^^*) $ dovc(^l,^^Co|^l,-B^,^o)| ^y , i=!/o uE , g i=Eoi - 

i 



(D13) 



In the formula for $dev, we have to suppress all informa- 
tion about the HF modes and that about response prop- 
erties of the LF mode. Preserved is only the information 
about the output current, 

$o(Co) = $dcv(0,0,0,0,Co|0,0,0,0,0) 

= (T:exp(<oJo):) . (D14) 

We also have to suppress the back-action of the detector 
on the amplifier and of the amplifier on the light source, 
which means dropping the corresponding differential op- 
erators, 

«{&£)Hi-£h* <D15) 

Consider firstly a photodetection formula in terms of the 
light source and composite detector. The corresponding 
part of Eq. (ID13|) reduces to, 

*o(Co)=$ s (^A^,-^A R o) 

x$ cd (( \E ,E*)\ Eo=0 , (D16) 

where 

<I> s fa,^o) = $dcvAfa,K)| : ) 

= (T:exp (iiSoVo - iv i>l) :) s , 
* cd (Co \E ,E*)= $ devSC (0, 0, 0, 0, Co | E , E* ,0,0,0) 
= (T: exp (i(oJo) :> cd - 

(D17) 

The differential operator in (|D16|) may be rewritten as a 
quantum average, 

^( r:eX K^0 AR02? ° + 4 A ^) : ) s 

= ( r:ex K^ 0+ A4); (D18) 



We used here the fact that, under the time- normal order- 
ing, classical radiation laws apply directly to operators, 
so that we could write, 

£ Q (t) = J dt'A m (t - t')V Q (t'). (D19) 

Similar to Eq. (|C20|) , 

= {T:T(£ ,£t):) s . (D20) 
Equation (|D16[) may therefore be written as, 

$ (Co) = (r:$cd(Co|^o,4 t )=) s - (° 21 ) 

Remembering the definitions of <£> and $ cd , this is equiv- 
alent to the two-layer average, 

(T:exp(iC ol 7 ):) o 

= (r:[(r:exp(zCoJo):) cd ]| £o£ ,^ o4 :), (D22) 

To decifer this relation, recall that averages denoted 
■ ■ • :) , are by definition conditional on Eo(t). In 
terms of current averages, 

(T:J (ti)---Jo(t m ):) o 
= <T: [(T:X(ti) • • • X(t m ):)J \ Eo ^ io ^, CD23) 

which is an under-the-RWA counterpart of Eq. ([75]) . 

The rest of Eq. (|D13[1 under approximations (|D14|1 . 
(|D15[) is manupulated similarly. So, photodetection 
statistics in terms of the properties of the three devices 
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reads, 



(T: exp (< X) :) G = <V: exp (^7^0 + 



x (T :exp («CoJo):) c 



(D24) 



E =Ei=0 



Again, recall that averages (T: • • • :) a and (T: ■ ■ ■ -)d are 
by definition conditional on, respectively, Eo(t) and 
Ei(t). Manipulations similar to those leading to Eq. 
(ID23I) allow one to write this relation as a three-layer 
average, 



<T:^ (tx) • • • ^o(^):> = <T:{<T= [(T:^^) - - • ^ (^):> d ] [ Ei> ^^ i; ^ =>J U ^,^ o ^ o t => s - (D26) 

I 



This relation fully justifies "doing quantum electrody- (T: ■ ■ ■ :) a and (7~: ■ ■ ■ -)d were subject to semiclassical 
namics while thinking classically" in Sec. IIVFI in the models, 
main body of the paper, where quantum averagings 
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